Two-component perfect fluid in FRW universe 



V.E. Kuzmichev, V.V. Kuzmichev 



Bogolyubov Institute for Theoretical Physics, 
National Academy of Sciences of Ukraine, Kiev, 03680 Ukraine 



(N 



O 



Abstract 



We consider the cosmological model which allows to describe on equal foot- 
ing the evolution of matter in the universe in the time interval from the inflation 
^ , till the domination of dark energy. The matter has a form of a two-component 

' perfect fluid imitated by homogeneous scalar fields between which there is en- 

£C) , ergy exchange. Dark energy is represented by the cosmological constant, which 

is supposed invariable during the whole evolution of the universe. The matter 
, changes its equation of state with time, so that the era of radiation domination 

O -1 in the early universe smoothly passes into the era of a pressureless gas, which 

^' ( , then passes into the late-time epoch, when the matter is represented by a gas of 

bj) 1 low-velocity cosmic strings. The inflationary phase is described as an analytic 

continuation of the energy density in the very early universe into the region 
^vq ' of small negative values of the parameter which characterizes typical time of 

, energy transfer from one matter component to another. The Hubble expansion 

CO ' rate, energy density of the matter, energy density parameter, and deceleration 

, parameter as functions of time are found. 

^ 1 PACS numbers: 98.80.Qc, 98.80.Cq, 95.35.+d, 95.36.+X 

: 

o ■ 

; 1 Introduction 

The standard cosmological model HQ El describes the universe which contains rela- 
tivistic matter (radiation), visible and dark matter, and dark energy in the form of 
the vacuum (cosmological constant) or quintessence (hypothetical cosmic fluid with 
the equation of state parameter which changes with time). The universe evolves so 
that in the course of the expansion, one or another matter-energy component begins 
to dominate over the others. 

In the early universe the dominating matter component is radiation with the 
energy density p ~ R~ 4 , where R is the cosmic scale factor. Then, as the uni- 
verse expands, the non-relativistic matter (dust) with the energy density p ~ R 
begins to dominate. If one assumes that in the late-time universe there exists a 
complementary linear constraint between the cosmic scale factor and total amount 
of matter in the universe as a whole, the latter phase will pass into the stage in 
which the matter in the universe will look like a perfect gas of low-velocity cosmic 
strings with p ~ R~ 2 (3j. In the radiation-dominated universe the number den- 
sity of photons is n 7 ~ -R~ 3 , and the energy of every photon decreases, during the 
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expansion of the universe, as m 7 ~ R due to the cosmological redshift. As a 
result, the effective mass of the matter attributed to relativistic matter reduces as 
well, M e jf ~ m^n^R? ~ R" 1 . Arguing in the same way, one finds that in the uni- 
verse, in which the matter is represented by a dust, the effective mass is constant, 
M e ff = const, expressing the constancy of the sum of the masses of the bodies in 
the volume ~ R 3 . In the low- velocity cosmic strings phase, the effective mass of 
the matter increases with the expansion of the universe, M e jf ~ R. On the time 
interval from the reheating after the inflation till the domination of dark energy, the 
equation of state of the dominating matter changes from the equation for radiation 
with the pressure p = ^pto the equation for the low- velocity cosmic strings p = — ^p 
passing through the pressureless state p = 0. 

It seems justified to consider the approach in which the equation of state of the 
matter in the universe changes continuously in time, passing through the limiting 
cases being specified. Models, in which continuous change of state of the matter 
was assumed, were studied by a number of authors [H[5]. 

In the present article we consider the model of the universe, in which the matter 
has a form of a two-component perfect fluid, while dark energy is represented by the 
cosmological constant, which is supposed invariable during the whole evolution of 
the universe. Components of a fluid are described by spatially homogeneous scalar 
fields. The first component has the equation of state which changes in time from 
the stiff Zel'dovich type equation to the vacuum type one. The second component 
is pressure-free matter. During the course of the evolution of the universe, there 
is redistribution of energy between the components. As a result, the total matter 
changes its equation of state in time, so that the era of radiation domination in 
the early universe smoothly passes into the era, when the matter is represented 
by a pressureless gas. The latter changes into the epoch in which the matter is 
described by a gas of low-velocity cosmic strings. The solution obtained for the era 
of radiation domination in the very early universe can be analytically continued into 
the region of extremely small negative values of the parameter which characterizes 
typical time of energy transfer from one matter component to another. At the same 
time, the first matter component acquires the properties of the vacuum with the 
corresponding equation of state and it can be identified with the inflaton field, used 
in models of inflation, in this stage. The contribution from the cosmological constant 
can be neglected in the early universe, while in the late-time epoch its role is vital 
to produce the accelerating expansion. 

For the homogeneous, isotropic and spatially flat universe with the Friedmann- 
Robertson- Walker (FRW) metric, in terms of a two-component model, we find the 
explicit expressions for the Hubble expansion rate, energy density of the matter, 
energy density parameter, and deceleration parameter as functions of time for the 
whole time interval. The limiting cases corresponding to small, intermediate (when 
the matter in the universe is represented by a dust), and large (the late-time epoch 
of low- velocity cosmic strings) values of time are considered for a comparison with 
known solutions. 
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2 Two- component perfect fluid 



Let us consider the FRW universe which is described by a set of equations 

tt2 ( 8ttG. 

* S [r) = — <" + 'A), (1) 
p + 3H(p + p) = 0, P a = 0, (2) 

p = W (t) p, p A = -p A , (3) 

where R(t) is the cosmic scale factor, pit) is the energy density of the matter which 
has a form of the homogeneous perfect fluid, p(t) is its pressure, w{t) is the equation 
of state parameter, p\ is the energy density of the vacuum which is connected with 
the cosmological constant A in the usual way as pa = g^jj, Pa is its pressure, G is 
the Newtonian gravitational constant, an overdot denotes d/dt, t is the proper time 
(units c = 1 are used). 

We assume that the matter consists of two components, 

P = Pb + Pd, P = Pb+ Pd- (4) 

It is supposed that a redistribution of energy between these components takes place 
in the course of the evolution of the universe. The energy conservation equation for 
the matter ([2]) rewritten for components takes the form 

p b + 2,H{p b +p b ) = Q, (>d + 3H(p d +p d ) = -Q, (5) 

where Q describes the interaction between the components. 

The components of the perfect fluid are imitated by scalar fields 4>b{t) and 4>d{t) 
with potentials V b {4>b) and Vd(cj>d), 

1-2 1-2 

Pi = ^4>i + Vi, Pi = -4>i -Vi, i = {b, d}. (6) 

The models of such a type which include a coupling between the matter components 
were considered in the literature, in particular, within the context of inflation and 
reheating and the coincidence problem of dark energy and matter in the acceler- 
ating universe (see, e.g., Refs. [SI [71 [8] and references therein). The form of the 
interaction term Q may be derived from different physical arguments or obtained as 
a solution of some dynamical equation, which describes the required properties of 
the matter fields fa. Since, according to Eqs. CO)-©, dark energy is modeled by the 
cosmological constant, only the evolution of the matter as a two-component perfect 
fluid is considered below. 

Let us assume that the field (pd forms the pressure- free matter component (dust), 

^4d=V d , p d = 2V d , Pd = 0. (7) 
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The potential Vd is taken equal to the density p&, Vd = Pb- Then, pd = 2pb, and the 
total energy density, the pressure, and the parameter w are 

p = 3p b , p = Pb, w = ^-. (8) 

In this case, Q = 2Hpb and the set of equations © reduces to one equation 

p b + 3H (p b + ^ Pl ) = 0. (9) 

The following condition is imposed on the second field (pb, 

^ b 2 = Vb e-^-^, (10) 

where to is the instant of time close to which the field (pb reproduces the pressure- free 
matter, r is the parameter which characterizes some typical time of energy transfer 
from one component to another and the value 1/r determines the mean rate of 
change of the equation of state parameter w(t) of Eq. ([5]). 

The field (pb, which satisfies the condition (|10p . describes the matter component 
which evolves in time from the stiff Zel'dovich matter, 

p b ps p b at t « 0, for — > 2, (11) 

r 

through the dust, 

p b w at t « t , (12) 
to the matter with vacuum-type equation of state, 

Pb ~ ~Pb at t > t - (13) 
Substituting Eq. (fTUj) into ([6]) and using the expression for w ([8]), we find 

w (t) = -Itanh^^^^) . (14) 

Thus, in the model under consideration, the equation of state parameter for the 
two-component matter w changes from ^, when Eq. (jlip holds, passing through the 
point w = for the pressure- free matter component (fl"2|) . to — | in the case (fl"3j) . 
taking all intermediate values. 

Using the mechanical analogy, the function (I14p can be considered as the antikink 
solution of the equation 

ltf + [-U(w)]=0, U = ^ (^-^) 2 > (15) 

which describes the motion of the analogue particle with zero energy in the potential 
[— U(w)] (cf., e.g., Ref. [9]). This potential has two maxima at the points w = ±| 
and a local minimum at w = 0. The analogue particle moves along the 'trajectory' 



4 



([HI) from the value w = 4 in the distant past (t = — oo) to the value w = — | 
reached at t = oo. At the moment t = to, the particle passes through the minimum 
of the potential at w = 0. Leaving the point w = 3, the analogue particle can only 
approach the point w = —3 at i — > 00, where its velocity and acceleration vanish. 
It cannot return back to w = i . This analogy with the motion of the particle allows 
to understand a unidirectional evolution of matter in the universe governed by Eq. 
([T]) from the radiation-dominated era through dust domination to the hypothetical 
late-time epoch, in which the matter is described by a gas of low-velocity cosmic 
strings. We note that Eq. (I15p has another solution in the form of the kink which 
is equal to the function (|14p with an inverse sign. This case would correspond to 
the model in which a gas of low- velocity cosmic strings at t = — 00 transforms into 
radiation at t = 00, but we do not consider it in the present paper. 

Let us find, how the cosmological parameters, namely the Hubble expansion rate 
and the total energy density, depend on time. From Eqs. ([1]), Q and (I14h . it follows 
the non-linear equation for the Hubble expansion rate, 



tanh 



- (t - 1 ) 

T 



9 A 

H 

3 



0. 



The general solution of this equation is 



H(t) 



coth Z (t) , 



where we denote 



Z(t) 



1 /A 



2 V 3 



Ct + 3t - rlncosh ( - (t - to) 



(16) 



(17) 



(18) 



C is the constant of integration. 

Substituting the solution (fT?]) into Eq. ([1]), we find the energy density of the 
matter as a function of time 



pit) = pa sinh 2 Z (t) . 
The energy density parameter of the matter, Qm = p/(p + pa), is 

Q M (t) = cosh' 2 Z (t) . 
The deceleration parameter, q = —1 — H/H 2 , is equal to 

3 - tanh [(t - t ) /r] 



-1+ 



(19) 
(20) 

(21) 



2 cosh 2 Z (t) 

The scalar fields and 4>d themselves can be restored by integrating the expres- 



sions 




Pit) 



Their potentials are 



e 2(t-t )/r + I ' 



- Idt^pjt). 



(22) 



(23) 
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3 Limiting cases 



Let us consider the different limiting cases, when the matter takes the forms of 
radiation (t <C t ), dust (t ~ to), and a perfect gas of low-velocity cosmic strings 
(t > t ). 

1. The late-time epoch of low-velocity cosmic strings. 

The Hubble expansion rate, the energy density, and the deceleration parameter 
in this case do not depend on the choice of the constant of integration C and the 
parameters to and r, 

H (t) = ^1 coth (JjtJ > = P A sinh_2 f/f *1 > 

H M (t) = cosh" 2 ' ?(*) = - tanh2 ^V^*) &t * ^ ^ 

If the cosmological constant is equal to zero, A = 0, then 

The equations PH) reproduce the result of Ref. [3]. The equation for p(t) (f25j) has 
a form of Whitrow- Randall's relation |10| . 

2. T/ze epoch, when the matter is in the form of dust. 
In this epoch the parameters are 

A . (\ [K ,_. „. \ . . _ 2 / 1 I A 



H{t) = y - coth ^-y - (3t + Ct ) J , p(t) = PA smh- 2 ^ - (3t + Ct ) j , 

Q M (t) = cosh- 2 ^iy|(3t + Ct )^ , ?(t) = -l + ~Q M (i) at t~t - (26) 

Choosing the constant of integration C = 0, one restores the expressions obtained 
in Ref. [31 [TTJ. In the case A = and C = 0, we have the following well-known 
result [12], 

3. T/te epoch of radiation. 
We have 

tf(t) = yfcoth [4t + (C - 1) t + Tin 2]^ , 

p(t) = PA sinh- 2 ^y?[4t + (C-l)t + rln2]j , 

M (i) = cosh- 2 ^y|[4t + (C-l)t +Tln2]j , 

g(i) = -l + 2Q M (i) att<Ci (28) 
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Choosing C = 1 — j- In 2, in the case of A = 0, we obtain the known result [12] 

In the special case, when t <C 4| (C — 1) to+ r In 2 1 (it includes the very early universe 
with t > M" 1 , where M p = G~ 1 / 2 (h = c = 1) is the Planck mass), it follows from 
Eq. (|28p that the energy density p remains almost constant and equals to 

Pe^Asinh- 2 ^y|[(C-l)io + rln2]J . (30) 

Passing to the limit A = 0, we obtain 

3 1 



2nGt 2 ( c _ 1+ t ln2 ^ 2 



to 



(31) 



4. T/ie era of inflation. 

Concerning the evolution of the matter, the era of inflation demands special 
examination. In the approach under consideration, we can perform an analytic 
continuation in the expression (|3ip into the region of small negative values of the 
parameter r, r < 0, |r| <C to- Then, it follows from Eq. (fTUj) that 4>b ps and 
the corresponding matter component has the vacuum-type equation of state, pb ~ 
—p b pa —Vb- It can be identified with the scalar field (inflaton) used in models of 
inflation to describe the evolution of the very early universe [21 [13]. Really, in this 
epoch, the energy density is p = 3pb ~ 3V;,, where 

Vb = 2^7 hi — ^ (32) 

plays the role of the inflaton potential. Assuming that the inflationary era is 
smoothly connected to the radiation-dominated phase, from Eqs. (f3T|) and ([32]) 
it follows that the constant C must be close to unity, so that the potential (|3"2"j) is 
very large. Setting C = 1 and |r| ~ Mr , we can rewrite Eq. (|32|) in a simple 
form V b w iM*. Substituting p w > p A into Eq. $1]), we find the solution 

8tt 



i?(t) ~ exp ^y^MptJ, which describes inflation. 

In the model under consideration, the inflationary era is driven by one of the 
components of the two-component matter, while the cosmological constant remains 
invariable at all stages of the evolution of the universe. In such an approach, the 
contribution from the cosmological constant can be neglected in the early universe, 
while in the late-time epoch of accelerating expansion its role becomes crucial, ac- 
cording to Eqs. ([24]) and (EST) . 



4 Discussion 

In Fig. 1, it is shown the time dependence of the decelerating parameter calculated 
in accordance with the general solution (|2ip and in cases when the relations for 
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Figure 1: The decelerating parameters (|2ip . (I24p . (j26[) . and (I28p versus time. The 

value of ^J~^ *o = 0.855 following from the WMAP7 data [14] is used. For definite- 
ness, the constants C = 1 and r = 0.2to are chosen. 

the epochs of low- velocity cosmic strings (|24p . dust (|26p . and radiation (|28p are 
extrapolated to the whole time interval. The free parameters A and to (as the age 
of the universe) were taken from the WMAP7 data [14] . All curves go to —1 for 
infinite t/to- The deceleration of the universe turns into acceleration at the instant 
of time t ~ 0.48to- If to is identified with the age of the universe, the corresponding 
cosmic redshift is z ~ 0.28. It can be seen from Fig. 1 that, near the instant of time 
t = to, the curves for a dust and radiation are close to the curve corresponding to 
the general solution (|2ip . whereas the curve for strings is far away from it. This 
can be interpreted as evidence that in the region t ~ to, effectively, the dominating 
matter component is a dust with a small addition of the radiation, while a gas of 
cosmic strings makes a negligible influence on the dynamics of the universe here. 

For a two-component perfect fluid with vanishing cosmological constant, A = 0, 
the solution of Eq. (|16p has a form 

(33) 

In this case the deceleration parameter does not depend on the constant of integra- 
tion C, 

g(t) = i{l-tanh[(t-t )/r]}. (34) 

This function is plotted in Fig 2. The deceleration parameter changes from 1 to as 
t/to changes from to infinity. Between the radiation-dominated phase (t/to 0.5) 



H{t) 



Ct + 3t - r In cosh ( - (t - t ) 
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Figure 2: The decelerating parameter (I34p versus time. It is taken r = 0.2to- 



and the era of cosmic string gas domination (t/to > 1-5), there is the transitional 
domain, where the phase of the radiation domination smoothly passes into the era 
of dust domination (t/to ~ 1)> which, in turn, then passes into the late-time epoch 
of cosmic string gas domination. Due to the gravitational attraction of the matter, 
the universe will be decelerating on the time interval from t/to = to t/to ~ 1-5. 
From the instant of time t/to ~ 1-5 to infinity, the effective gravitational attraction 
of the perfect gas of low-velocity cosmic strings will be practically compensated by 
its negative pressure acting repulsively, and such a universe will expand with the 
constant velocity. In the region of the large values of t/to, the cosmic scale factor 
depends on time as R ~ t. As it is claimed, e.g. in Ref. [15j . power-law cosmology, 
R ~ t a with a ~ 1, is consistent with the recent Wilkinson Microwave Anisotropy 
Probe (WMAP7), Baryon Acoustic Oscillations (BAO) and Hubble data, therefore 
the consideration of models, in which the cosmic scale factor evolves linearly with 
time, while the cosmological constant is negligible, may be of practical interest. 



Comparing Fig. 1 with Fig. 2, one finds that the general behavior of the decelera- 
tion parameter with time remains the same, however, the action of the cosmological 
constant smoothes the curve and also produce a new physical effect which results in 
a change of the sign of the deceleration parameter, so that the deceleration phase in 
the history of the universe passes into subsequent period of accelerating expansion. 
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5 Conclusion 



In the present paper, we propose a single approach to the description of the evolution 
of the matter in the universe with non-zero cosmological constant in the whole time 
interval from the inflation till the domination of dark energy. The new element here 
is the inclusion of the hypothetical late-time epoch, when the matter is represented 
by a gas of low-velocity cosmic strings, in the general scheme. In some sense, these 
strings may be considered as emergent ones, since they correspond to the collective 
motion of the matter. We draw attention to the fact that, if the final stage of 
the evolution of the matter in the universe is not a pressure-free matter, but a gas 
of low-velocity cosmic strings, then there arises an additional 'symmetry' in the 
equation of state of the matter. Namely, the equation of state parameter w (|14p 
changes from ^, passing through the point w = 0, to — |, taking all intermediate 
values. For such a model, we find the solution (|17p of the equation for the Hubble 
expansion rate (|16D in analytical form as a function of time. The Hubble expansion 
rate (fT7|) , the energy density of the matter (fT9j) , and the deceleration parameter ([2T]) 
as functions of time reproduce all known cosmological solutions (f25l) , (|27|) , and ([29]) 
for the different limiting cases obtained by a number of authors [3j \10\ [IT\ [T2] for 
vanishing cosmological constant. 
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